Soft nuclear equation-of-state from heavy-ion data and implications for compact stars 
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Measurements of kaon production at subthreshold energies in heavy-ion collisions point to a soft 
nuclear equation-of-state for densities up to 2-3 times nuclear matter saturation density. We apply 
these results to study the implications on compact star properties, especially in the context of 
the recent measurement of the two solar mass pulsar PSR J1614-2230. The implications are two- 
fold: Firstly, the heavy-ion results constrain nuclear matter at densities relevant to light neutron 
stars. Hence, a radius measurement could provide information about the density dependence of 
the symmetry energy which is a crucial quantity in nuclear physics. Secondly, the information 
on the nucleon potential obtained from the analysis of the heavy-ion data can be combined with 
restrictions from causality on the nuclear equation-of-state. From this we can derive a limit for the 
highest allowed compact star mass of three solar masses. 

PACS numbers: 26.60.-c Nuclear matter aspects of neutron stars 26.60.Kp EoS of neutron stars 



I. INTRODUCTION 

Recently, the measurement of the Shapiro delay for 
the millisecond pulsar PSR J1614-2230 provided a new 
reliable limit for the highest known pulsar mass of 
(1.97 ± 0.04) M [T]. Though, also higher neutron star 
masses up to (2.74 ± 0.2) M Q for PSR J1748-2021 are 
under discussion [5J [5] , these are not fully obtained from 
the observation of effects from general relativity and 
are therefore less reliable. The value for the maximum 
neutron star mass is of major interest for nuclear physics 
as it is tightly connected to the properties of dense 
nuclear matter, especially the stiffness of the nuclear 
matter equation-of-state (EoS). 

The exploration of nuclear matter at high density 
and temperature is also a topic of large experimental 
effort. High energy nucleus-nucleus collisions at collider 
facilities like RHIC at Brookhaven and the LHC at 
CERN explore the nature of strongly interacting matter 
at large temperatures, while the FAIR facility at GSI, 
Darmstadt, and the NICA facility at Dubna will focus 
on matter at high baryon densities. 

Today, firm information exists about the properties 
of nuclear matter around its saturation density uq. 
The compression modulus Kq = 9^ |„ describes 
the stiffness of isospin symmetric nuclear matter at 
n Q ~ 0.17fm~ 3 . It has been determined to (235±14)MeV 
from the study of giant monopole resonance oscillations 
of heavy nuclei [4]. The second important quantity 
is the symmetry energy S(n). Ground state masses 
of nuclei, the neutron skin thickness of heavy nuclei, 
and measurements of giant dipole resonances point to 
S(n ) = So ~ (28-32)MeV 0|6]. Isospin diffusion data 
from intermediate energy heavy-ion collisions provides a 



constraint on the density dependence of the symmetry 
energy around n ~ n a 7 . The properties of nuclear 
matter beyond no can be extracted from collective flow 
of nucleons in non-central nucleus-nucleus collisions [5]. 
However, by selecting non-central collisions, the probed 
density enhancement is less pronounced. Furthermore, 
the flow analysis with respect to the compression mod- 
ulus of nuclear matter is not yet conclusive [5] and the 
results obtained by different transport models do not 
agree [5]. 

The focus of this work lies on the impact of a soft 
nuclear matter EoS, as obtained from the KaoS exper- 
iment at GSI, Darmstadt, for baryon densities up to 
(2 — 3) n n on the properties of compact stars. We make 
use of results on K+ meson production in nuclear colli- 
sions at subthreshold energies, wherein the probability to 
produce K + is highly dependent on the reached particle 
density in the collision zone, most pronounced at the 
lowest bombarding energies. The level of compression 
is in turn controlled by the stiffness of nuclear matter 
through the nucleon potential Uat. The more attractive 
Uat is, the higher is the produced K + meson abundance, 
which therefore can serve as a probe for the stiffness of 
nuclear matter |10j . 

This paper is organized as follows. First, we review 
the results on K + meson production from the KaoS 
experiment and the need of a soft EoS to explain 
the observed data. In the next section, we study the 
relevance of the KaoS results for low mass neutron stars 
and their connection to the nuclear symmetry energy. 
This question arises inevitably since light neutron 
stars have low interior densities which can be in the 
same range as the ones probed by KaoS. As will be 
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discussed in the following, the knowledge of the isospin 
symmetric nuclear EoS from the KaoS experiment could 
enable a direct study of the symmetry energy from light 
neutron stars. The last part of the paper is focused 
on the determination of the upper limit on neutron 
star maximum masses via the highest possible mass. 
This can be deduced applying the stiffest causal EoS 
starting from a fiducial density which is set by the 
reliability of the nuclear EoS for densities below the 
critical value. As the highest possible neutron star mass 
is anti-proportional to the critical density, a large value 
of the latter provides stronger restrictions. Previous 
publications discuss values for the highest possible mass 
assuming critical densities beyond hq. However, as we 
will argue, our study is the first one to use a nuclear EoS 
which is experimentally obtained at supra-saturation 
densities and thereby allows the choice of a fiducial 
density larger than the saturation density. We conclude 
the paper with a discussion of the obtained limits for 
the highest possible neutron star mass with respect to 
neutron star observations, such as the PSR J1748-2021 
pulsar. 



n ~ (2 — 3) no [TH Q21 [TS1 [TB], the KaoS measurements 
probe nuclear matter in this density range [TH [17] . As 
a consequence, the stiffness parameter K < 200 MeV 
characterizes the nuclear EoS for n ~ (2 — 3) no while the 
compression modulus Kq — (235 ±14) MeV is determined 
by the properties of nuclear matter at no- 
The measured attractive nucleon potential and the cor- 
responding soft nuclear EoS [TTHT3] have to be tested 
on their compatibility with neutron star observations. 
Establishing a link between nuclear properties obtained 
from experiments and neutron star observations is a non- 
trivial task [I5H2D] . Nevertheless, in this work we want to 
investigate the implications of the results from the KaoS 
collaboration on compact stars. To avoid any assump- 
tions about the unknown high density nuclear matter, 
we focus on neutron star properties which can be directly 
linked to the soft EoS at n ~ (2 — 3)no. The first chosen 
observables are the radii and moments of inertia of light 
neutron stars, which will be discussed in the following. 
In the last section of the paper, we focus on the limit for 
the highest allowed compact star mass [21] . 



II. A SOFT NUCLEAR EQUATION-OF-STATE 
FROM SUBTHRESHOLD K+ MESON 
PRODUCTION IN A+A COLLISIONS 

The measurements on K + meson production in nuclear 
collisions at subthreshold energies were performed with 
the Kaon Spectrometer (KaoS) at GSI. The beam energy 
dependence of the K + multiplicity ratio from Au+Au and 
C+C collisions (M/A)a u +Au/(M/A)c+c at subthreshold 
energies of 0.8 to 1.5GeV per nucleon has been introduced 
as a sensitive and robust probe for the stiffness of nuclear 
matter [IT] . 

To describe the experimental results, IQMD and RQMD 
(Isospin and Relativistic Quantum Molecular Dynamics) 
transport model calculations were performed [121 113] . 
applying a Skyrme type Un with two- and three-body 
forces. Two parameter sets were chosen so as to repro- 
duce two nucleon potentials with different levels of repul- 
sion at supra-saturation densities. The transport calcula- 
tions consistently demonstrate that the beam energy de- 
pendence of the kaon multiplicity ratio is described best 
by the nucleon potential with small repulsion. When ap- 
plied to infinite isospin symmetric nuclear matter the lat- 
ter corresponds to a soft nuclear equation of state charac- 
terized by a stiffness parameter K = 9^ |„ = 200 MeV 
[T2"] . While the KaoS measurements seem to suggest even 
lower values of K < 200 MeV 12 , it should also be noted 
that the isospin asymmetry of the colliding Au+Au sys- 
tem could result in an increase of K by ~ 10% [T4] . 
Despite being defined in the same way, the stiffness pa- 
rameter K and the compression modulus Ko must be 
clearly distinguished from each other. Since the K + 
mesons in the Au+Au system are found to be pro- 
duced predominantly at supra-saturation densities of 



III. IMPLICATIONS OF THE NUCLEAR 
SYMMETRY ENERGY ON RADII AND 
MOMENTS OF INERTIA OF LOW MASS 
NEUTRON STARS 

With their immense interest in nuclear physics [22j . 
neutron star radii are a focus of major attention |23H30j . 
Limits on radii are obtained from X-ray bursters and low- 
mass X-ray binaries (see [24] [25] I30H36] and references 
therein) whereas such studies could be substantially re- 
fined in the future with the Large Observatory For X-ray 
Timing (LOFT) [37]. Recently, Bauswein and Janka [28] 
have studied the tight correlation between the frequency 
peak of the postmerger gravitational-wave emission with 
radius measurements, while Guillemot et al. |38j have 
presented a new technique to deduce limits on the mo- 
ment of inertia of gamma-ray pulsars. 
Neutron star radii and moments of inertia are tightly 
connected to the symmetry energy |39j . very similar to 
the nuclear symmetry energy dependence of the neutron 
skin thickness of heavy nuclei [3D]. For our study light 
neutron stars are of special interest due to their low cen- 
tral densities. If the latter are in the same range as the 
one probed by the KaoS experiment, information on the 
stiffness of nuclear matter in the neutron star interior can 
be taken directly from the KaoS data. Being the remain- 
ing uncertainty of the nuclear EoS, the symmetry energy 
can then be extracted by radius and moment of inertia 
measurements of light neutron stars. 
For a consistent analysis, neutron star properties and the 
KaoS data should both be described by the same under- 
lying EoS. We choose a Skyrme type EoS with a nucleon 
potential similar to the one which was used in the anal- 
ysis of the KaoS data. The energy per baryon is written 
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as 
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(u§ - F(u)J + S> 7 ] , (1) 



whereas u = n/no is the baryon number density and 
Eq is the average kinetic energy of nuclear matter at uq 
and a proton fraction of Y p = 0.5. Two- and three-body 
forces, described by the terms B and D, together with 
cr, are fitted to reproduce the binding energy per baryon 
E(riQ,Y p = 0.5) = —16 MeV, the stiffness parameter 
K, and the saturation density uq. While So is varied 
between 28 MeV and 32 MeV, its density dependence 
is chosen as a power law with u 7 where 7 = 0.5 — 1.1. 
The values for Sq and 7 are motivated by heavy-ion 
experiments and recent analysis of neutron stars in 
X-ray binaries El EH 02] . Eq. Q is used to describe 
matter in the neutron star core, while an inner and outer 
crust are incorporated following HI] . 
We first study neutron stars with M = 1.25 M Q , in 
accord with the lightest pulsar masses deduced so far 
from observations [45] . The radii R and moments 
of inertia I are calculated for non-rotating neutron 
stars, whereas their dependence on K is shown in 
Fig. [T] for different symmetry energy configurations. It 
can be seen that both largely depend on the density 
dependence of the symmetry energy in form of 7, i.e. 
L = 3n {dS(n b )/dn b )\n . At if ~ 200 MeV, stiff and 
soft symmetry energy configurations lead to a difference 
in the neutron star radius and moment of inertia of up 
to AR ~ 1.5 km and AI ~ 2.5 • 10 43 g cm 2 , respectively. 
The central densities of the corresponding stars are in 
the range of < 3.4 n , similar to the density region 
explored by the KaoS collaboration. At these densities, 
hyperons, kaons, and quarks, if they appear at all, 
should not play a dominant role [46j . However, their 
possible appearance in a first order phase transition 
could be probed in the neutrino signal of a galactic 
supernova explosion [47l [48] . 

In Fig. [2] we study the radii and moments of inertia 
of light neutron stars with masses in the range of 
(1.1 - 1.6) M Q for K ~ 200 MeV and with equal 
configuration for the symmetry energy as in Fig. [T] 
It can be seen that neutron stars with masses up to 
M < 1.3 Mq have central densities which are low enough 
for the entire neutron star interior to be described by an 
EoS which is probed by the KaoS experiment. With the 
isospin symmetric part of the nuclear EoS determined by 
the latter, radius and moment of inertia measurements 
from such stars could thereby distinguish between a soft 
and a stiff behavior of the symmetry energy. 
A high precision in neutron star radius measurements 
could be achieved for compact stars in accreting binaries 
by future observatories such as the proposed LOFT 
)37j . The latter has the promising potential to resolve 
neutron star radii with a 5% accuracy and would thereby 
be able to distinguish between a stiff and soft density 
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FIG. 1: Radii R and moments of inertia / of non-rotating 
neutron star configurations with 1.25 Mq in dependence of K 
for different setups of So and 7 from Eq. (1), i.e. different 
slopes of the symmetry energy L. Central densities are given 
in units of the nuclear matter saturation density no. 



dependence of the symmetry energy. For neutron star 
moments of inertia, only few suggestions exist to deduce 
informations from stars in specific systems |49j . Guille- 
mot et al. 3* suggest a new method to obtain a lower 
limit on / of gamma-ray pulsars using the spin-down of 
the star coupled to the gamma-ray efficiency. Though 
this approach requires a very good knowledge of the 
star's distance, it seems to be especially interesting for 
low mass neutron stars. Due to the specific shape of 
the exclusion curve the latter appears to provide larger 
restrictions to stars with a low moment of inertia. 

It is noteworthy that for this study we use a sim- 
ple ansatz for the nuclear symmetry energy which is 
usually taken to study neutron stars and the density de- 
pendence of the symmetry energy in heavy ion collisions 
[5]. To make precise predictions of light neutron star 
radii and moments of inertia it is of course necessary to 
know the exact form of the nuclear symmetry energy, 
e.g. also a possible non-monotonic behavior with density 
[50] . In addition, large rotational frequencies can 
increase the neutron star radius and moment of inertia, 
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FIG. 2: Radii and moments of inertia of light neutron stars 
for Ko ~ 200 MeV in dependence of their mass for different 
symmetry energy setups as in Fig. [I] Central densities are 
given in units of the nuclear matter saturation density no. 



whereas the latter can also be affected by superfluidity 
in the neutron star interior. 

Nevertheless, our study demonstrates that a direct 
cross-check between results from heavy ion experiments, 
such as KaoS, and astrophysical observations, in form 
of radii and/or moments of inertia of light neutron 
stars, is possible, since both probe nuclear matter in 
the same density region. Moreover, with the isospin 
symmetric part of the nuclear EoS at supra-saturation 
determined from the KaoS data, we provide predictions 
concerning the relation between the compactness of light 
neutron stars and the density dependence of the nuclear 
symmetry energy. 



IV. MAXIMUM NEUTRON STAR MASSES 

We turn now from low mass stars to the maximally 
allowed gravitational mass of a compact star. The 
most massive stable neutron star configuration depends 
strongly on the stiffness of nuclear matter. The softer 
the nuclear EoS, the lower is the maximum mass which 



can be reached before the star becomes unstable against 
collapse to a black hole. 

Ref. [H] introduced the idea to use known properties of 
hadronic matter up to a fiducial density ricru- At higher 
densities, the low density EoS is smoothly connected to 
the stiffest possible EoS allowed by causality. It is given 
by 



p(e) =Pcrit + (e- e C rit), 



(2) 



where p is the pressure and e the energy density |51j . 
Their values e cr u and p cr it are taken at the transition 
density n cr it from the known low density EoS to the 
stiffest one. The causal EoS exerts the highest pressure 
to balance the gravity of the neutron star mass and there- 
fore leads to the highest possible maximum mass. The 
existence of such a stiff EoS beyond (2 — 3) no can be as- 
sociated with the emergence of new physics, such as the 
appearance of quark matter [55], and could be probed 
by the CBM experiment at FAIR in the near future [53] . 
Previous studies showed that the highest possible mass 
can be expressed as [35J [50] [5U |55J : 



M high =4.1M (e crit /e ) 



-1/2 



(3) 



Hereby, the applied EoSs for hadronic matter were ob- 
tained from measurements of nuclei at n ~ n and then 
extrapolated to higher densities [54Tl56| . We want to 
point out that for such an approach the reliable fiducial 
density can only be taken as 



Ecrit — £0 



n ~ 1.4 • 10 14 g/cm 



(4) 



As can be seen from Eq. ([3]), this results in Mhigh ~ 
4.1 M . Since for higher e crit the nuclear EoS is deter- 
mined via extrapolation from eo, the corresponding lower 
values of M hig h, e.g. M hig h ~ 2.9 M for 

^crit — 2 6q [54 

[56] . represent only possible predictions for the highest 
possible mass, and depend on the validity of the extrapo- 
lation. Results from the KaoS collaboration, on the other 
hand, directly test the EoS for densities n ~ (2 — 3) no 
and therefore justify the choice of n cr i t in the same range. 
Consequently, in this section we will calculate the highest 
possible masses with the restriction that for n cr it up to 
(2 — 3) no the nucleon potential should fulfill the results 
from the analysis of the KaoS data. For Skyrme type 
models this restriction corresponds to K < 200 MeV. 
Today, there exists a large variety of Skyrme parame- 
ter sets [57]. For our work we chose two representative 
approaches for Skyrme models, BSk8 and SLy4 [551 159j . 
which are applied in astrophysical and nuclear physics 
studies. Both models are fitted to reproduce properties of 
neutron rich nuclei and nuclear matter at saturation den- 
sity. They have similar values of K ~ 230 MeV while the 
symmetry energies have a different density dependence. 
For Relativistic Mean Field (RMF) EoSs, the stiffness 
of high density nuclear matter is determined by the nu- 
cleon effective mass m* at no and not by K [5D] . Conse- 
quently, the corresponding Un is chosen by varying m* 
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Therefore, we find that while the BSk8, SLy4, and NL3 
parametrizations are fitted to nuclear matter at satura- 
tion density they are not applicable for higher values of 
n ~ (2 — 3)hq. However, for comparison, we keep these 
parameter sets in the following calculation. 
Fig. [4] shows the highest possible neutron star masses for 
n crit ~ (2— 3)no and the discussed EoSs. For n crit ~ 2hq, 
the star is dominated by the stiffest EoS and therefore 
reaches masses of up to 3 M Q . Smaller maximum masses 
of Mhigh ~ 2.4 M Q are obtained for the upper limit of 
n crit ~ 3 n . The higher n cr i t is, the later is the onset 
of the stiffest EoS in the star's interior. Consequently, 
less mass is supported and the value for the highest mass 
decreases. 

It can be seen from Fig. [4] that for the Skyrme type EoSs 
(Skyrme K = 200 MeV, BSk8 or SLy4) as well as the 
BHF calculations the maximum masses are smaller than 
for the TM1 or NL3 models. While having a nucleon 
potential well below the KaoS limit, RMF models with 
varying m* give a maximum mass which is above the 
ones for the Skyrme and BHF EoSs. The smaller m* 
is, the stiffer is the EoS and, therefore, the larger is the 
maximum mass. We checked the impact of symmetry en- 
ergy on our results, varying So and 7 for Eq.([T]) and the 
p meson coupling constant in the RMF models. Since 
the maximum mass configuration is dominated by the 
causal high density EoS and not by asymmetry at low 
densities, the symmetry energy has very little influence 
on Mhigh, leading to a difference in mass of the order 
AM < 0.02 M Q . As can be seen from Fig. |4j we can 
conclude that a pulsar mass of 2.7 M Q as found by [2] 
is marginally compatible with a soft EoS and requires a 
prompt transition from a soft EoS to the stiffest possible 
at a density around (2.2 — 2.5) n . 



for a given Kq, so as to obtain a nucleon potential which 
is similar to or even more attractive than the Skyrme 
parametrization within the density limits. A more at- 
tractive nucleon potential at supra-saturation densities 
allows a higher compression of matter for the same bom- 
barding energy, which enhances multiple-step processes 
in subthreshold kaon production and is in line with KaoS 
results. 

Fig. [3] shows the nucleon potentials of the chosen mod- 
els together with the restriction from the KaoS experi- 
ment. Realistic parameter sets should have a Um below 
the KaoS constraint up to densities of n ~ (2 — 3) no- 
The TM1 {61J parametrization as well as the Brueckner- 
Hartree-Fock approximation (BHF) 62-65 fulfill this re- 
quirement. For a RMF model with no = 0.17 fm -3 and 
K = 220 MeV, we arrive at m*/m = (0.53 - 0.65) 
for n ~ (2 — 3) uq while for Ko = 250 MeV we obtain 
?7i* jra = (0.54 — 0.67) for the same density range. Other 
schemes, such as BSk8 [58 , SLy4 [59] or NL3 [66], pro- 
duce nucleon potentials which are more repulsive than 
the Skyrme benchmark in the density region of interest. 



V. SUMMARY AND CONCLUSIONS 

K + multiplicities from heavy-ion collisions at GSI indi- 
cate that the nuclear EoS is soft for densities of (2 — 3) no- 
To test the consequences on neutron star properties we 
study low mass neutron stars whose radii and moments 
of inertia are sensitive to the density dependence of the 
symmetry energy. We find that light neutron stars with 
M < 1.25 Mq are strong candidates to have central den- 
sities in the range of ~ 3 710 which is within the density 
region tested by the KaoS collaboration. Such low mass 
neutron stars are therefore well suited for future radius 
measurements. 

Moreover, to test whether a soft nuclear EoS up to ~ 3 no 
is compatible with tentatively massive neutron stars such 
as PSR J1748-2021 [2], we apply the KaoS results up 
to densities of 2 n < n crit < 3 Uq |21j and introduce 
the stiffest possible EoS for n > n cr i t to calculate the 
highest allowed maximum mass Mhigh- Within this ap- 
proach, the KaoS results confirm the previous theoretical 
estimation for the highest possible neutron star mass of 
~ 3M Q . A pulsar mass of 2.7 M Q is not excluded by the 
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KaoS data, but requires the onset of the stiffest possible 
EoS already at n cr i t ~ 2.2 — 2.5 no- In the future, the 
CBM experiment at FAIR will probe densities beyond 
3 no using rare probes as e.g. D-meson [53] and thereby 
provide further constraints on Mhigh- 
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